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Abstract: 

We construct explicit models of particle dark matter where the attractive force in the 
dark matter sector creates a narrow near-threshold resonance that qualitatively changes 
the energy dependence of the annihilation cross section. In these models, the resonant 
enhancement of the dark matter annihilation can easily source the excess of energetic 
leptons observed by experiments on PAMELA and FERMI satellites. The distinct feature 
of these models is that by construction the enhancement of the annihilation cross section 
shuts off when the dark matter velocity falls below the typical Milky Way values, thus 
automatically satisfying constraints on dark matter annihilation imposed by the CMB 
anisotropics and gamma ray constraints from satellite galaxies. 
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1. 


Introduction 




It 


is well-known that baryonic matter contributes only about 2( 


)% of the mass density in the 



Universe. A plethora of independent cosmological observations strongly support that the 
remaining 80% is composed of cold dark matter (DM) . The most recent observation suggests 
that DM contributes a fraction of Odm — 0.229 ± 0.015 to the total energy content of the 
universe [1] . Notwithstanding our rather precise knowledge of the global dark matter energy 
density, its nature remains a mystery, and constitutes one of the most profound questions 
in fundamental physics. Several candidates have been proposed to account for DM [2], with 
weakly-interacting massive particle (WIMP) featured prominently among them, thanks to 
the natural mechanism for the its relic abundance regulated through thermal annihilation 
cross section of about (av) ~ 3 x 10~ 26 cm 2 s _1 . The annihilation of DM into the Standard 
Model (SM) particles in the local Universe, and specifically inside our galaxy, can result 
in its "indirect detection" via energetic annihilation products such as high energy gamma 
rays, charged leptons, protons and and anti-protons. 

Recent measurements of cosmic ray spectra of charged leptons by the PAMELA, ATIC, 
Fermi and H.E.S.S experiments have reported some "anomalies" relative to the prior ex- 
pectations based on theoretical models of cosmic rays. In particular, the PAMELA results 
show a rise in the relative contribution of positron flux into the total e + e~ flux above 
energies of around 10 GeV to 100 GeV [3]. Also, a broad excess in the total e + e~ spectrum 
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at around TeV scale energies was reported by several experiments [4, 5, 6]. On the other 
hand, the flux of antiprotons does not show any anomaly relative to prior expectations [7]. 
While it is possible that some unaccounted astrophysical sources are behind the excess of 
energetic leptons, the DM explanation to these excesses has also being entertained, steering 
theoretical modifications of the minimal WIMP paradigm. These excesses are not easily 
accommodated in the thermal WIMP scenario, because of the need to enhance the annihi- 
lation cross section into leptons in the galactic environment by a large factor O(10 2 — 10 3 ), 
while maintaing the antiproton fraction within limits. Both model-building goals, enhance- 
ment of the cross section, and suppression of antiprotons is naturally achieved in models 
that supplies the DM with the the so called " dark force" , a relatively light compared to the 
WIMP mass mediator particle that connects DM with the SM [8, 9, 10, 11]. The enhance- 
ment of the annihilation rate in galactic environment can be broadly described as due to 
the extra mediator-exchange attractive force that creates the Coulomb (or Sommerfeld) 
enhancement of the cross sections at low WIMP velocities. Ref. [9] summarizes certain 
broad classes of enhancement mechanisms within this scenario: 1. the usual Sommerfeld 
type enhancement that can be approximated as Seoul — ^a' ' jv in the limit of small mass 
of the mediator, with v being WIMP relative velocity and a' the coupling constant of 
the dark force; 2. narrow resonances that can bring the enhancement far above 5couij 3. 
recombination into WIMP-anti-WIMP bound state with the emission of a mediator, that 
brings further numerical enhancement to the cross section. 

Further developments included an in-depth calculations of the Sommerfeld enhancement 
as the function of the mass of the mediator and possible small mass splittings of WIMP 
states [12], and a discussion on whether these models afford the required enhancement 
factors [13, 14, 15, 16]. It is fair to say at this point that simple versions of such scenarios 
are under strong pressure from precision data on anisotropics of the cosmic microwave 
background (CMB) from WMAP7 [17, 18], and the (non) observation of gamma rays from 
dwarf spheroidal galaxies by Fermi-LAT [19, 20, 21, 22, 23]. Neither of these latter studies 
have so far given any indication on an enhanced DM annihilation signal. However, the 
kinetic energy of WIMP particles in dwarf spheroidals and in the universe during the CMB 
decoupling is much smaller than the average kinetic energy of DM in our galaxy. The 
momentum of DM redshifts in the same rate as the expansion of the Universe after its 
thermal decoupling. Therefore, the average kinetic energy of DM at the recombination 
must be smaller than the temperature of the universe, which is about 0.3 eV. Inside dwarf 
spheroidal galaxies, the velocity dispersion is around or smaller than 10 km s _1 [24, 25], 
which is one order of magnitude smaller than 200 km s _1 , typical for DM velocities in our 
galaxy in the ~ kpc vicinity of the Solar System. Therefore, if the enhancement factor 
depends on the kinetic energy of dark matter in a non-monotonic way and peaks when the 
relative velocity of DM is around 200 km s _1 , it can co-exist with constraints from the 
CMB and dwarf spheroidal galaxies, and yet generate the desired enhancement aimed at 
explaining the anomalous spectra of positrons and electrons with DM annihilation. 

The purpose of our work is to show that the existence of narrow near-threshold resonances 
in WIMP- WIMP system can easily invalidate the extrapolation of galactic annihilation rate 
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Figure 1: Schematic illustration of the main idea about the importance of narrow near-threshold 
resonance for the DM phenomenology. In black, are the av in units of pb for a generic Sommerfeld- 
enhanced annihilation with the saturation at low energies and the narrow p-wave resonant cross 
section with the peak at E = 1 MeV. The annihilation rate is given by the integral from the overlap 
of av with the DM energy distribution, which are pictured as Maxwell distributions for energetic 
particles in the early universe (blue), galactic environment (red) and sub-haloes (green). Narrow 
resonance can have a large overlap with galactic DM distribution, while minimizing the rate both 
for large and small velocity dispersions. Sommerfeld enhanced model is guaranteed to have the 
enhancement for the small velocity dispersions. 

to lower velocities, enabling to escape constraints from the CMB and dwarf spheroidals. 
The main idea of our work is illustrated pictorially in Fig. 1. The narrow resonance (black 
curve pointing upward at 1 MeV, T ~ 0.1 keV) maximizes the overlap with DM kinetic 
energy distribution inside the galaxy, represented on the plot by the Maxwell distribution 
with T = 1 MeV. At the same time, the overlap with the early universe distribution (blue 
curve), and less energetic distribution for dwarf spheroidal galaxies (green cure, T = 0.01 
MeV) is minimized. This should be contrasted with the case of a typical Sommerfeld- 
enhanced cross section (black curve with flat behaviour at E — ¥ 0, Fig. 1) that has 
significant overlap with both red and green distributions. 

A general formula for the resonance-enhanced scattering cross section can be written 

as 



where k is the momentum in the centre-of-mass frame, J is the total angular momentum 
of the resonance, S\ and S2 are the spins of two incoming particles, E is the centre-of-mass 
energy, Tin and r out are incoming and outgoing widths. Notice that since the resonance 
is close to the zero-energy threshold, the entrance width must depend on the momentum 
of the particles and can be parameterized as r; n (E) = r{°\E/E R ) l+1 / 2 , where I is the 
angular momentum of the system. In the limit of the very narrow resonance the energy 
dependence of Eq. (1.1) is often approximated by the delta- function. 
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The origin of the resonances may vary. It can be generated by the intermediate s-channel 
mediator with the mass closely tuned to 2m\yiMP> Ref. [26]. A more natural, i.e. less tuned, 
scenario might be that two DM particles form a resonance through a long range attractive 
force [27, 9, 28]. While the standard Sommerfeld enhancement models do exhibit a series 
of resonances (see e.g. [12]), these are not suitable for our purposes because of their large 
width. Therefore, a dedicated exercise is required in order to find whether scenario in Fig. 
1 can have concrete model-building realizations. In this paper, we propose two prototype 
models, which can induce an O(10 3 ) resonant enhancement factor for the DM annihilation 
in our local galaxy, while avoiding constraints from the CMB and dwarf spheroidal galaxies. 
We view these models as a "proof of existence": not overly complicated models of dark 
sector can indeed furnish the required non-monotonic behaviour of the enhancement factor, 
which calls for more in-depth analysis of phenomenology of such models. 

The rest of the paper will be organized as follows: In Sec. 2 and Sec. 3 we propose two 
candidate models that have narrow resonances either in s or p wave. In Sec. 4 we provide 
a discussion of the gamma-ray constraints from the galactic centre. We summarize the 
results in Sec. 5. 

2. Narrow resonance in the s-wave 
2.1 The model 

The main idea of the model is to realize the s-wave resonance via the small mass splitting 
of two WIMP mass eigenstates, xi aR d X2, with X2 being the lightest and stable. The 
following annihilation chain is envisaged: 

X2 + X2 -> (xiX2)i« resonance ->■ light mediators -)• SM. (2.1) 

While our models can be constructed for different spin of dark matter particles, we prefer 
to work with scalars as the most economical and simplest option. 

We choose the Lagrangian of the dark sector to consist of new £7(1) gauge group, one 
complex DM scalar x an d three sets of Higgs fields, H, A and A c : 

C = ^^(A ff /ff 2 +Lc.)-i(A A AV+Lc.)-i(A^t x 2 +h.c.)- J(A AAc A 2 A c2 +h.c.) , 

(2.2) 

where Cq includes the SM Lagrangian, the new £7 (1) field Lagrangian as well as other terms 
that are not crucial for the DM phenomenology. Under the dark £7(1) symmetry x an d H 
are singly charged, while A and A c cary the charge ±2. We assume that this dark £7 (1) 
symmetry is spontaneously broken by the vacuum expectation value (vev) of A c . Without 
loss of generality, we assume coupling constants A aa c , A A and Xh to be real and leave only 
Aa to be complex. Then after spontaneous symmetry breaking, the real and imaginary 
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parts of x and A are split into the mass eigenstates, with eigenvalues given by 
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The residual interaction Lagrangian between x and A can be written as 

1 

1 



V2 



Ai 
A 2 



^7T^(Xi - X2) ~ ImA A XiX2 



ImA A 2 2n . t, x 
— ^— (Xi - X2) + ReA A XiX2 



(2.4) 



In the primordial Universe, the thermal annihilation of dark matter stops when the temper- 
ature drops to a level T stop ~ M x /20. The mass difference between the real and imaginary 
parts of x is negligible compared to T stop , and one can treat xi and X2 as one complex 
scalar. To discuss the thermal annihilation, let us consider the case with ReA# ~ ImXn ^> 
other four-point couplings. Then, there are two main annihilation channels, H scalars and 
vector bosons, with the following cross sections: 
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(2.5) 



and its sum will determine the DM abundance. For the choice of parameters considered 
below, the H channel is the dominant one. 

2.2 Formation of the bound state 



The ladder exchange by leads to the attractive force between xi and X2 particles, Fig. 
2. It also enables the attraction in the X1X1 X2X2 channel. One can readily observe 
that U-exchange respects certain charge symmetry in XiXj sector: z + j-even and z + j-odd 
sectors do not mix until one includes Higgs-induced interactions (2.4). The attractive V- 
mediated force leads to the bound states in both sectors, but by increasing the mass of V or 
reducing the couplings, the bound states can be squeezed above X2X2 separation threshold. 
We choose the parameters in such a way that X1X2 state is just above the threshold. 

Let us consider the first chain in process (2.1), namely the capture into (xiX2)ls resonance. 
From d x t\ one can see that both Ai and A2 induce a process of X2X2 - > XiX2- Since in the 
absence of the Higgs-mediated interaction the two sectors, X2X2 and X1X2 do not mix, it is 
rather easy to obtain a very small entrance width for this resonance. After some lengthy 
but straightforward calculations, we derive 
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where ob = 27r/ (ayM x ) is the Bohr radius of the bound state and v r is the relative velocity 
between the two \2 particles in the continuum. 

We next calculated r out in our model. It turns out that X1X2 pair annihilation into one 
gauge boson and one scalar in the final state is negligible. (It is either in p-wave of is 
additionally suppressed by a small vev). Since the annihilation happens at short distances, 
one can conveniently tie r ou t to the annihilation cross section of X1X2 — > H scalars: 



( av )xiX2->HH = {<JV) 

so that the out-going width becomes 

Tout ~ (™)o4tt|V(0)| 2 = 
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(2.8) 



It is important to specify the ratio of two widths in our model: 



r in _ 2v r ReA A ImA A 
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(2.9) 



In order to have a large enhancement and a narrow width at the same time one would like 
to have 

r in ~ T out ~ keV , (2.10) 

which removes a lot of freedom in the choice of parameters. Substituting explicit expressions 
for qb into r out , we have 



Tout ~ -a^M^Ho » afy GeV 



(2.11) 



For the WIMP mass of 1 TeV, the desired value for r ou t is achieved with ay ~ 0.01, so 
that ayM x ~ 10 GeV and a# ~ (5 GeV)" 1 . Moreover, the requirement of the primordial 
annihilation rate to be consistent with the present day DM energy density, (<tv)q ~ 1 pb, 
further imposes the requirement of A# ~ 1. Therefore, to make Ti n « r out we need 
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(2.12) 



Since Ai and A2 both mediate an attractive interaction for the xi pair, we have to make 
sure that this interaction is not generating strong Sommerfeld-type enhancement so that 
the model cannot be ruled out by satellite galaxies and CMB constraints. In the region 
m A (1 2 ) S> 5 GeV and 3> m\ 2 the above relation can be rewritten as 




(2.13) 
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Figure 2: Diagrams illustrates the formation of bound state. 



If one defines a dimensionless coupling in the following way 

the Sommerfeld enhancement factor generated by A2 can be written as 

s ^ _ 7ra A M x _ A A A a ^ 0.3A A (2 15) 

rrtA 2 2m a 2 16M X M x 

Therefore, one can readily see that as long as Aa < My, A2-exchange will not generate 
any enhancement in the annihilation of two X2 particles. 

We now address the Sommerfeld enhancement due to exchange of V, as shown in Fig. 2(b). 
The effective potential for |xiXi) and \X2X2) system can be written as [8, 12] 





The potential is attractive for one of the combination of IxiXi) and \x2X2) and repulsive for 
the other. However, in the case that 2(M Xl —M X2 ) > -©binding) where ©binding is the binding 
energy of the Is bound state, one cannot get a sizeable Sommerfeld enhancement. 

Picking up a specific point in the parameter space with M X2 — M Xl = 8 MeV (eg = 
x/2(M X2 - M Xl )/M x /a v » 0.4), ©binding = 7 MeV, and m v = 2.6 GeV, in Fig. 3 we show 
the dark matter annihilation cross section as a function of the kinetic energy E^, including 
both X2X2 and X1X2 intermediate channels. The blue dotted curve is the contribution from 
the X1X2 resonance, while the green dashed curve is the contribution from the Sommerfeld 
enhanced X2X2 annihilation. 

Assuming the distribution over kinetic energy to be Maxwellian, we also plot (av) as a 
function of temperature, Fig. 4. One can readily observe the desired behaviour anticipated 
in the Introduction: while the cross sections are strongly enhanced for an O(MeV) temper- 
atures, they are abruptly diminished to the primordial values as soon as the temperature 
drops below 50 keV. 

Finally we discuss the relation between my and ay that places X1X2 resonance just above 
the X2X2 separation threshold. It can be derived from the Schrodinger equation 

V 2 

y + V(r)$> = E$>, (2.17) 

2fM x 
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Figure 3: av as a function of kinetic energy of dark matter. The blue dotted and the green dashed 
curves show the contributions from the X1X2 resonance and the X1X1 annihilation 




0.001 0.005 0.010 0.050 0.100 0.500 1.000 

T (MeV) 



Figure 4: (av) as a function of temperature. The red, green and blue curves are for E r = 0.5, 1, 1.5 
MeV, respectively. One can observe a sharp drop-off in the averaged annihilation cross section below 
T = 50 keV for all three cases. 



where V(r) = —^ e - m vr an( j a _ |_ anc j ^ ~ M x /2 is the reduced mass of the X1X2 
system, and E < since it is a bound state (for the discussion of my — ay relation, the 
transition to X2X2 continuum is a small perturbation and can be neglected). 

From the straightforward numerical analysis of this equation we derive the relation between 
my/(ayM x ) and E/(a v my), shown in Fig. 5. The resulting domain of parameters is 
consistent with the requirement my/(ayM x ) < 0.6 that we can obtain from our analysis 
01 r in (out)- 
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Figure 5: Relation between my/(ayM x ) and \E\/(ctyM x ). 



3. Narrow resonance in the p-wave 
3.1 The model 

In this model we assume that the dark matter particle is a complex scalar interacting with 
two gauge bosons and the effective Lagrangian can be written as 



We would like V\ to generate a 2p resonance during the annihilation of dark matter when 
the kinetic energy is around 1 MeV, and we would like the mass of V 2 to be smaller than the 
energy difference between the 2p resonance and the Is bound state so that the 2p resonance 
can decay to Is state. However, in order for avoid overly large enhancement of the cross 
section due to light V 2 exchange constrained by the CMB, we need g 2 coupling to be small. 
Thus, the main idea behind this model is to realize the p-wave resonance-mediated capture 
into the sub-threshold Is state: 

X + X f -> (xX t )2 P resonance -> (xX^is bound state + V 2 ->• light mediators ->• SM. (3.2) 
In this model, in the early universe the main process for dark matter annihilation is 



(3.1) 



(3-3) 



with the cross section given by 



av = 




(3.4) 



where notation a.\ = 4<7^/(47r) is introduced. The relation between M x and a.\ follows from 
the normalization to standard WIMP annihilation cross section. 
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Figure 6: The blue solid and red dashed curves show the numerical calculation and T[ n = aE 3 / 2 
fit for sin 2 S as a function of energy, for the choice of M x = 1000 GeV, ot\ = 0.025, myi = 3 GeV. 



3.2 2p bound state, Tin and T 
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At low energies the dynamics of a WIMP anti-WIMP pair is governed by the attractive 



Yukawa potential V(r) 



-m vi r 



jr. We analyze the Schrodinger equation in this 



potential numerically imposing the condition of near-threshold 2p resonance. Our study 
shows that for M x = 1 TeV and ay = 0.025 the choice of ray j '(avM x ) = 0.1201 places 
this resonance at E = 1 MeV. One can also quantify the degree of fine tuning required 
for having an MeV resonance by normalizing it on the binding energy of the 2p state in 
hydrogen-like case, Er/(o{ M x /16) ~ (1 — few)%. 

We then perform the study of the elastic scattering on this resonance to determine Tin. 
The cross section for p-wave resonance scattering can be written as 
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with the explicit dependence of Tin on energy, r; n = aE 3 ^ 2 . The coefficient of proportion- 
ality a can be calculated numerically, and for the same choice of parameters (M x = 1 TeV, 
ot\ = 0.025, my i = 3 GeV) we determine a = 1201 GeV -1 / 2 . The resulting phase shift as 
a function of energy is plotted in Fig. 6. 

We next address the out-going width that has two contributions, r ou t = F 2p annihilation + 
r2p-s>is+V2 : one from the annihilation of x an d directly in p-wave, and the other from 
2p — > Is decay with emission of V%. The Is WIMP-onium state of x an d X* immedi- 
ately decays, but its rate of decay is inconsequential for the size of the annihilation cross 
section. 

For the reference, the the width of Is bound state is given by 

T ls « (at;)o47r|V>(0)| 2 « 1.25a? GeV « 20 keV . (3.6) 
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Figure 7: av around the p-wave resonance. The red, green and blue curves are for r out = 10, 1, 
0.1 keV, respectively. 

Since the direct annihilation in a 2p resonance is further suppressed compared to the above 
width by a factor of 0(af), we conclude that 

T2p annihilation < 10 eV , (3.7) 

which is much smaller than what we need. According to the discussions in the previous 
section, in order to generate enough enhancement, r ou t has to be around 1 keV. Therefore, 
decay width from the 2p resonance to the Is bound state should be around 1 keV, which 
is easy to achieve via 

Tout ~ r 2p ->i s +y 2 - f-J a.2a\M x . (3.8) 

After fixing the T ou t, we can use Eq. (1.1) to calculate av around the resonance, which is 
shown in Fig. 7. 

From Eq. (1.1), we can further infer that in our case r ou t <C T[ n , and the total width is 
quite narrow at the resonance position, so that (av) can be written as 

r 

l x u o 

where vq ~ 10~ 3 is the average velocity of DM in the local galaxy. Eq. (3.9) suggests that 
in order to get enough enhancement r ou t should scale as 0.1 x (M x /1 TeV) 3 keV. From 
the relic abundance relation (3.4) we can further infer that 

n/lTeV\ 2 

a 2 ~ 1.3 x 10^ 2 — — . (3.10 
V M x J 

Assuming again Maxwellian distribution, the averaged av is shown in Fig. 8. As the one 
can see, besides the p-w&ve induced resonance, the cross section also has a "flat part" 
associated with the Sommerfeld enhancement of the s-channel annihilation to be discussed 
below. 



(av) res ~ — ^ , (3.9) 
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Figure 8: Thermally averaged av as a function of "temperature" assuming the velocity distribution 
is Maxwellian. The red, green and blue curves are for M x = 3000, 5000, 7000 GeV with r out = 0.5, 
3.3 and 10 keV, respectively. 

3.3 s-channel enhancement and constraints from CMB 

The vector boson V2 may induce a large s-wave enhancement for the annihilation between 
X and when the kinetic energy is small, which will be strongly constrained by the 
observation of CMB. At the era of last scattering, the kinetic energy of the dark matter is 
much smaller than eV so that the dependence of the Sommerfeld enhancement factor on the 
kinetic energy of dark matter must be negligible [8] . Fig. 9 shows the s-wave enhancement 
factor induced by V\ and V2 as a function of M x in the condition that a± is fixed by the 
thermal relic abundance and 02 is determined by r out from Eq. (3.10). 

We quote recent interpretation of WMAP7 data as a limit on the dark matter annihilation 
cross section at the 95% C.L. [17]: 



where / is the energy deposition factor. Assumption of / ~ 1 means that most of the energy 
from the annihilation of final states is efficiently transferred to the primordial plasma. 
Taking / = 1, we show the CMB sensitivity by the blue line in Fig. 9. We can see that in 
the region with M x above 4 TeV the model is generally safe from the CMB constraints, i.e. 
the 2p resonance can enhance local annihilation, while the enhancement in the s-channel 
is within limits. 

4. Constraint from annihilation in galactic center 



/ 



(<">)cmb 2.42 x 1(T 27 cm 3 /s 
M x < GeV 



(3.11) 



The models presented in the two previous sections avoid the CMB and dwarf spheroidal 
gamma constraints by design. This, however, does not exhaust all constraints on dark 
matter annihilation, and in this section we address constraints that usually come from 
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Figure 9: The s-wave Sommerfeld enhancement factor induced by V\ and V2 as a function of M x 
with a.y\ fixed by the requiring that the process XX ~^ V\Vi generates the right relic abundance of 
dark matter in the early universe by thermal annihilation and a\/2 is determined from Eq. (3.10). 

(non) observation of gamma excess from the central region of our galaxy. For concreteness, 
we shall follow recent works that discuss implications of H.E.S.S. data on high energy 
gamma-rays from the galactic centre (GC) [29, 30]. One familiar difficulty is that the 
density profile of dark matter in the GC is quite uncertain. In our model of dark resonance, 
this uncertainty is compounded by the velocity dependence of the cross section {crv), which 
brings the sensitivity of the rates to the distribution of dark matter over velocities in GC. 
Given this velocity dependence of {crv), we argue that there is no conflict between these 
constraints and the enhancement factor for the DM annihilation within ~1 kpc from the 
Solar System, relevant for the PAMELA signal explanation. 

The H.E.S.S. telescope is sensitive to high energy gamma-rays from a few hundred GeV 
to a few tens of TeV, and can use its superior angular resolution to select the gamma-ray 
signal from the GC, and compare its signal with the "background region" . The background 
region is chosen to be farther from the GC than the signal region in order to capitalize 
on the enhancement of the DM density in the signal region relative to the background 
region, inferred from numerical simulations of the DM profiles [31, 32, 33]. Consequently, 
the average number of gamma-rays from the signal region should be larger than from the 
control region, and the difference between the two can be used for setting upper limits on 
the annihilation cross section. 

This technique is very sensitive to the assumptions about the density profiles, as emphasized 
in Ref. [29]. For example, a constant-density core of DM within the inner region of around 
450 pc of the Milky Way will wash away all constraints that can be derived from HESS 
GC analysis. We now investigate how the velocity dependence of (cry) will modify the 
conclusions of analysis in Refs. [29, 30]. For concreteness, we shall assume a NFW density 
distribution for the Milky Way. 
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Figure 10: Left: NFW distribution, the solid and dashed curves show the ratios of velocity 
dispersion and density distribution normalized to the same quantity near the Solar System. The 
two shaded regions are the signal region and control region, respectively. Right: the distribution of 



(crv) in the signal region and control region in the second model with M x — 5000 GeV, r o 
keV. The solid and dashed parts show the signal region and control region, respectively. 



3.3 



The differential flux per solid angle generated by the DM annihilation can be written 

as 
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where p is the energy density of DM, r ga i = y Rq — 2xRq cos(Z) cos(6) + x 2 is the distance 
between the line-of-sight point x and GC, Jo is the normalization constant defined as 
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The difference of our analysis from Refs. [29, 30] is that here (av) also depends on r ga i 
through the dependence on the velocity distribution. 



The NFW profile is given by the expression 

p(r) = 



Po 



(1 + 



(4.3) 



Ra 



where for the Milky Way R s is about 20 kpc. Then the velocity dispersion can be calcu- 
lated from Jeans equation and the numerical solutions is shown in Fig. 10. More details 
on the velocity structure of the central region can be found in [34]. For simplicity, we 
further assume that the velocity distribution is Maxwellian. Then the velocity dispersion, 
normalized on the dispersion around the Solar System, shown in Fig. 10, can be rigidly 
linked to the average of kinetic energy, and can be further interpreted in terms of effective 
temperature. Therefore, from Fig. 10, we can see that in the signal region the effective 
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temperature is reduced by about one order of magnitude whereas only by a factor of a few 
in the control region. Therefore, in models with the resonance "optimized" for the anni- 
hilation around the Solar System, the ratio J^ ff (signal) / (control) becomes smaller and 
the ensuing constraint is weaker. Since the information given in Ref. [29] is not sufficient 
for an independent analysis, we only checked pixel and pixel 2 in Fig. 2 of Ref. [29], 
which are in the signal region and control region, respectively. Without considering the 
velocity dependence of (<jv) we get 



which are close to the average values in the signal region and control regions. Then, after 
including the velocity dependence of the annihilation cross section, we can get 



Therefore, we can see that after considering the velocity dependence, the contrast between 
control and signal region is eroded, and the constraints from the H.E.S.S. GC gamma-ray 
observation become much weaker. 

5. Summary and discussions 

Models with enhanced annihilation rate of dark matter aimed at explaining the PAMELA 
positron excess and Fermi-LAT e + e~ spectrum are strongly constrained by the absence 
of any evidence for non-thermal energy injection in the CMB data, and by the absence of 
gamma-ray excess from dwarf spheroidal galaxies and the galactic centre by Fermi-LAT and 
H.E.S.S. respectively. In particular, models that have a Sommerfeld-type enhancement of 
annihilation at small velocities, that scales as inverse velocity within some range, generally 
predict that the annihilation cross section (cju)cmb, dw.sph. > (f"w)gaiaxy I n this paper, we 
have shown that this is not the case in models where the annihilation proceeds via a narrow 
resonance at kinetic energy around 1 MeV. Indeed, the narrowness of the resonance allows 
to maximize the cross section right in the desired range of WIMP kinetic energies. 

We propose two prototype models. In the first one, the mass splitting in the WIMP 
sector was used to create the s-resonance. The narrowness of the resonance, crucial for 
the success of this construction, comes from the fact that two sectors, (X2X2) and (X1X2) 
exist independently, and do not mix unless one introduces a new Higgs interaction that 
connect these two sectors. Consequently, the width of X1X2 resonance can tuned to be 
around 1 keV, which automatically shuts down the enhancement of av when the kinetic 
energy of the DM particles becomes much smaller than an MeV. In the second model, we 
considered the p-wave resonance in the WIMP-aniti-WIMP state, so that the incoming 
width is proportion to where v is the velocity of DM. As a result, the low- velocity 
annihilation rate in this channel is automatically suppressed. It is fair to mention that 
these models are tuned, as the position of the resonance above the threshold is typically 
less than the energy level spacing in the bound state system, but we find that the degree 



•/pixel = 1015Jo , Jpixel 2 = 613Jq , 



(4.4) 



jpixel = 161 Jo , Jpixel 2 = 147Jq . 



(4.5) 
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of fine-tuning is 0(1%), and far less if one were to achieve a similar resonance with the 
exchange by a heavy particle in the s-channel. 

The two models found here do not exhaust all possibilities, as the mechanisms for achieving 
the narrow resonance in the two- WIMP state are quite numerous. Another generic mech- 
anism that we would like to mention here without building an explicit model for it goes 
as follows: the resonance can be achieved with the combination of attractive and repulsive 
force. For example, a Higgs and a vector exchange will create the following potential for the 
XX system: V xx = (ay exp(— myr) — ah exp(— mh,r))/r, while both pieces in the potential 
will be attractive for XX^ ■ With my < and ay < a^, there will be a potential minimum 
at small r, separated from large r by the barrier. Needless to say that this is exactly the 
case that can give narrow resonances, and the annihilation of XX state via such resonance 
can be achieved via a doubly-charged Higgs state. 

Given that narrow resonances in the two- WIMP system can qualitatively change the out- 
comes of the DM annihilation constraints, we hope that the full analyses based on the 
formula (1.1) with appropriate energy dependence of r; n will be performed in the near 
future. We close our paper with a series of additional comments: 

• In our scenario, the annihilation cross section of DM peaked at MeV scale, which 
is about the starting temperature for the big bang nucleosynthesis (BBN). Thus 
one expect the enhancement of the annihilation rate during the BBN, with possible 
observable outcomes, depending on the specifics of the final state particles at the end 
of the annihilation cascade. We leave this subject for the future study [35]. 

• Diffuse galactic gamma ray constraint will provide a strong limit on the models 
discussed here, because it is less sensitive to the distribution over the velocity. The 
most recent paper by FERMI collaboration [36] that appeared after our work has 
been completed, has shown to disfavour large enhancement of the annihilation rates 
for the TeV scale WIMPs. A few hundred GeV WIMP could still be entertained as 
the PAMELA signal explanation, and we believe that with some modifications, our 
models with the near-threshold resonance can be made consistent with 0(100) GeV 
WIMP mass scale. 

• While the enhancement of the lepton yield via a resonant channel was the primary 
reason for our investigation, we note that very similar resonant mechanisms can be 
invoked for enhancing the annihilation to monochromatic photons. The models of 
this type were previously discussed in [27]. 

• Models discussed in this paper do have a considerable degree of self-interaction, 
which can be constrained by the ellipticity of DM halos [13, 37], and the dynam- 
ics of the so-called "bullet" cluster [38]. The constraints from ellipticity is about 
two orders of magnitude stronger than from bullet cluster, and is given by jj- < 
4.4 x 10 -27 cm 2 GeV -1 . For TeV dark matter, it is equivalent to the sensitivity to 
elastic cross section at ~ 10 -24 cm 2 , which is considerably larger than the unitarity 
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bound in I = 0, 1 channels, thus not making any additional restrictions on models 
discussed here. 
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